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Abstract. For classical canonical transformations, one can, using the Wigner transfor- 
r^ . mation, pass from their representation in Hilbert space to a kernel in phase space. In this 

Oh| paper it will be discussed how the time-dependence of the uncertainties of the corresponding 

time-dependent quantum problems can be incorporated into this formalism. 
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^ ! 1 Introduction 



In classical Hamiltonian mechanics the time-evolution of a physical system is described by canon- 
ical transformations in phase space that keep the Poisson brackets of the transformed coordi- 



f^ ■ nate and momentum with respect to the initial ones unchanged. This transformation in phase 

OO . space can be described (for a one-dimensional problem in physical space and, therefore, a two- 

^■"■^ ' dimensional one in phase space, to which we will restrict ourselves in the following) by the 

so-called two-dimensional real symplectic group 5*^(2, R), represented by 2 x 2 matrices with 
^ ■ determinant equal to 1. (In order to compare our new results with earlier ones for the time- 

H \ independent case, we only consider the homogeneous symplectic group without translations, not 

the inhomogeneous symplectic group /S'p(2,IR).) It has been shown in [1] how it is possible to 
obtain the representation of the group of linear canonical transformations in time-independent 
quantum mechanics via the determination of the configuration space representation of the uni- 
tary operator that connects quantum mechanically the transformed variables x and p with the 
initial ones, x' and p' . A subsequent Wigner transformation shows explicitly that for the time- 
independent problems considered by this method, essentially the classical results are reproduced. 
This agrees with the fact that at least for quadratic Hamiltonians, the Wigner function evolves 
as W{x' ,p' ,t) = W{xm{x' ,p' , —t),pMix',p', —t),0), where xm and pM are the Moyal time evo- 
lution of position and momentum which, again for quadratic Hamiltonians, coincide with the 
classical evolution [21 [3l HJ [5] . 

In a different study [6] it has been shown that, for time-dependent quantum systems, cha- 
racteristic differences compared with the classical situation arise, in particular, when the time- 
dependence of the quantum mechanical uncertainties of position and momentum are taken into 
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account. One such difference can already be found for the most simple physical system, the free 
motion. 

In classical mechanics, the dynamics of a system may not only be described by its trajectory 
but, particularly for systems with many degrees of freedom, a statistical description in terms of 
virtual ensembles or distribution functions in phase space (also called F-space in this context) is 
possible. These density- type functions qy{x,p, t) must be globally conserved, which is guaranteed 
if they fulfil a kind of conservation law in the form of a continuity equation in phase space 
that connects the explicit change in time of the density function ^r with the divergence of a 
current jy = Qrvr- If the system obeys Hamilton's equations of motion, the divergence of the 
velocity field vr always vanishes, Vrvr = (where Vr is the nabla-operator in phase space). 
In a hydro dynamical description, this would mean that ^r describes an incompressible medium. 
In this probabilistic context, this corresponds to Liouville's theorem that, from all possible 
transformations of phase space, selects only those where a phase extension always retains its 
volume during motion. Making use of the concepts of measure theory, this statement can be 
made even more precise in the formulation that the measure of point-sets is an invariant of the 
time-evolution of the virtual ensemble. 

It can be shown straightforwardly that the quantum mechanical density ^(x,t)=^*(x,t)^(x,t), 
corresponding to the complex solution ^(x,^) of the time-dependent Schrodinger equation also 
fulfils a continuity equation (now in position- or configuration space). In this case, however, 
the divergence of the corresponding vector field is proportional to the relative change in time 
of the mean square deviation of position (position uncertainty) (x^) = (x^) — {x)'^, i.e., VxV = 
^^(a;^)/(2;^). So, VxV only vanishes if (x^) is constant which, e.g. for a Gaussian wave packet, 
would correspond to a constant width of this wave packet. It is, however, well known that 
already for the free motion, the wave packet width is not constant but spreading with increasing 
time. If one would consider the wave packet width as a kind of measure for the "volume" of our 
quantum system, it certainly would not be conserved under time-evolution. This situation also 
does not change if one tries to include the momentum aspect by also taking into account the 
corresponding mean-square deviation (p^) and then considering the product of (x^) and (p^) as 
an adequate measure for the "volume" of the system; since (p^) is constant for the free motion, 
therefore, the product with (x^) still grows in time. (For a more detailed discussion, see also [6].) 
A more consistent approach for a comparison with the classical case could, however, make use 
of the phase-space formulation of quantum mechanics in the form of the time-dependent Wigner 
function which shall be done in this paper. 

Unlike in the time-independent situation where the quantum mechanical results mainly re- 
produce the classical ones, in the time-dependent case there are obviously, at least, formal 
differences between the classical and the quantum mechanical descriptions of the system even 
already for such simple ones like the free motion. These differences are intimately connected 
with the time-dependence of the typical quantum mechanical aspect of the system, namely, 
the uncertainties of position and momentum. Therefore, in this paper, we will investigate the 
influence of the time-dependence of the quantum system, in particular of the uncertainties, on 
the representation of the group of linear canonical transformations in quantum mechanics. 

For this purpose, in Section^ we briefly summarize the main results of the time-independent 
quantum mechanical case. In Section [3] we then consider the time-dependent case and discuss 
the characteristic differences compared with the time-independent situation, in particular the 
role of the time-dependence of the quantum uncertainties. In order to be on the safe side 
of systems with exact analytic solutions, we restrict our discussion to systems with at most 
quadratic Hamiltonians, in particular, to the harmonic oscillator (with possibly time-dependent 
frequency) and the free motion. Possible ways of overcoming these limitations and of extending 
our method to further problems will be mentioned in Section [J] where the results will also be 
summarized and some perspectives will be mentioned. 
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2 Time-independent case 

The time-evohition in classical Hamiltonian mechanics is described by canonical transformations 
in phase space that can be represented by 

X \ fa h \ f x' 



where a, b, c, d are real and the determinant of the 2x2 matrix is 1, i.e., ad — be = 1. The 
group of transformations represented by the 2x2 matrices is the so-called two-dimensional real 
symplectic group 5p(2,M). 

Following [T] (see Chapter 35 and references cited therein), it is possible to obtain the repre- 
sentation of the group of linear canonical transformations ([1]) in quantum mechanics. Referring 
to [HI?], the main objective is to determine the configuration space representation 

{x\U\x') =K{x,x') 

of the unitary operator U that provides the quantum mechanical relation between x, p and x' , p' , 
according to 

x = Ux'U-^, p = UpU-^. (2) 

With the help of the kernel K(x,x'), the effect of any canonical transformation ([T|) can be 
described as 

r+oo 

*(x) = / dx' K{x,x')'^{x'). 



The integral kernel K{x, x') has been derived taking into consideration the fact that it must 
satisfy the following two differential equations [U [9] 



ax + b——- I K(x, x) = x'K(x, x'), 
I ox J 



(3) 



cx + d-^) K{x,x') = -^l^K{x,x'). (4) 

I ox J I ox' 

An exponential ansatz, bilinear in x and x', finally leads to K{x,x') in the form 

/Mi 



^^^'^'^ " \2^b) «^P |-^[«^' - 2xx' + dx'2] I . (5) 

This kernel K{x,x'), related with the specific canonical transformation, is formulated in 
configuration space, whereas, the corresponding classical canonical transformation is formulated 
in phase space. Therefore, it is interesting to discuss the representation of this canonical trans- 
formation in the phase space version of quantum mechanics that was developed by Wigner [10] 
with the help of the corresponding distribution function. This distribution function W{x^p) can 
be obtained from a given wave function ^{x) in configuration space via the so-called Wigner 
transformation 

"'<^-''' = 2^ /r •" "'"""'' (^ + f) * (^ - f) ■ "" 

Applying this transformation to the kernel (0) leads to the phase space kernel in the form [7] 

K{x^ x' ,p,p') = 5[x' — {ax + bp)]6[p' — {ex + dp)], 

showing that, for this linear canonical transformation, the kernel coincides with its classical 
limit. So, the quantum mechanical problem mainly reproduces the classical situation without 
any additional specific quantum mechanical aspect. 
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3 Time-dependent case 

Now we will investigate how far this is still true in the case of specific quantum dynamical 
aspects entering the problem. So, considering time-dependent problems in quantum mechanics 
in terms of the time-dependent Schrodinger equation or equivalent formulations, one finds that 
not only classical position and momentum change in time (in a way that can be described by 
canonical transformations) but, also the typical quantum mechanical degrees of freedom, like 
position- and momentum- uncertainties, may be time-dependent (corresponding, e.g., to wave 
packets with time-dependent width). For certain problems with exact analytical solutions in 
form of Gaussian wave packets, it has been shown (see, e.g., [11] and references cited therein) 
how the transition from initial position and time (in configuration space) to any later position 
and time can be achieved with the help of a time-dependent kernel (or propagator) K{x, x' , t, t') 
according to 

r+oo 

^{x,t)= dx' K{x,x\t,t')^{x',t'). (7) 



The integral kernel K{x,x' ,t,t') can be obtained in different ways, e.g., via Feynman's path 
integral method [12], or, for kernels quadratic in x and x' , similar to the discussion in the 
time-independent case [7] in the form 



^^"'"''*'*' = °)=i2;^Ito^J ^"Pl2M 



1 



/ m, \ 2 I im 



zx^ — 2x { — ] + u 



OLQ \ao 



/ \ 2 
X 



(8) 



where a^ = \J '"^ ° is proportional to the initial position uncertainty, or initial mean square 
deviation in space, (a;^)o = (3:^)0 — (x)q (where (• • • ) denotes mean values calculated according to 
(• • • ) = J_ °° \I/* • • • $ dx'). We choose this form of the kernel because of the explicit appearance 
of ao and, therefore, the initial position uncertainty. This quantity is essential for the time- 
evolution of the quantum uncertainties. As has been shown in [11], already for the simple 
example of the harmonic oscillator, the time-evolution of the position uncertainty and, hence, 
the width of the Gaussian wave packet solution of this problem, behaves qualitatively very 
different if the initial position uncertainty corresponds to that of the ground state (constant 
width) or differs from it (oscillating width) where the latter represents the general solution 
which, in the limit u; ^ 0, provides the free motion wave packet, whereas, the former only leads 
to a plane wave solution in this limit (for further details see [H]). The time-dependence enters 
this kernel explicitly via the parameters z{t) and u{t). In the limit f — > 0, the kernel turns into 
a delta function. 

Since, according to ([7]), the dependence of \I'(a;, t) on x and t enters only via if (x, x', t, t'), this 
kernel also must fulfil the time-dependent Schrodinger equation. Inserting K, in the form given 
in dS]), into this equation shows that the parameters z{t) and u{t) both fulfil the equation of 
motion for the corresponding classical problem (e.g., the free motion or the harmonic oscillator 
with possibly time-dependent frequency), however, they are not independent of each other but 
are coupled via the relation 

zu — uz = 1. (9) 

This resembles the condition that the determinant of the entries of the 2x2 matrix of the 
linear canonical transformation ([1]) must fulfil. 

Inserting ^ into ([7]) with an initial Gaussian wave packet ^y^p{x' ,t') of the form 



, , , , / mf3n \ * I im 



x'^ 


)%2^x' 


aoy 


/ m 
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leads to a Gaussian wave packet 

(10) 



^WP(x,t) = (-) U^ exp - 



z z{u + iz) 



whose maximum follows the classical trajectory 7?(t), which is, therefore, up to a constant, 
proportional to z, i.e., 

-(x)(t) = — -r?(i) (11) 



"oPo aoPo 

with initial momentum po and mean value of position (x) which is, according to Ehrenfest's 
theorem, identical with the classical trajectory ??(t). 

The coefficient of the term quadratic in x in the exponent of the Gaussian can be expressed 
via the complex quantity 

2h k 1 X 

m z z[u + iz) A 

where u and z have been combined to form the complex variable 

\ = u + iz. (13) 

The quantity —y fulfils the complex nonlinear Riccati equation (here given for the harmonic 
oscillator) 

—y+(—y) +uj^ = o, 14 

which can be linearized, using equation p2|) . to a complex Newtonian equation 

\ + uj'^\ = Q (15) 

for the variable A (for details see \\.1\ I13j). The complex quantity \{t) can also be expressed in 
polar coordinates as 

\ = ae^"^ = acosif + iasmip, (16) 

where relation ([9]) turns into 

^ = \, (17) 

which has similarities with the conservation of angular momentum but, here, for the motion 
of A(t) in the complex plane. 

Through equations (fT3]) , (fTSl) and (fTTI) it is obvious that the imaginary part of equation (fTSJ) , 
up to a constant factor, is just the Newtonian equation of motion for the wave packet maximum. 
The width of the wave packet (jlOp . or its position uncertainty (x^) = (x^) — (x)^, respectively, is 
connected with the imaginary part of (^y) via "^yi = 2m(x'^) ' -^'^o™ ^ i^ polar coordinates (see 
equation (fT6|) ) and pursuing equation (fT2]) and pT|) it follows that the variable of the complex 
Riccati equation (J14p can be written as 

2h \ \ a d 1 f2h \ f2h 

— y ] = T = - + w = - + ^^2 = \ ~yR + M ~yi 

m J A a a a^ \m J \m 
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with o? = 2m{Sp') /h. The imaginary part of equation (|14p agrees with equations (J17p and ([9]). 
The real part of equation p4|) with the above-mentioned results leads, however, to a real non- 
linear differential equation for the position uncertainty a{t), 

a + u?it)a = ^r. (18) 

It is well-known in the literature [14^ [TS] that a pair of differential equations consisting of 
the Newtonian equation of motion (jlSp . but now only for r/(t), and the nonlinear equation (jlSp . 
possesses a dynamical invariant, the so-called Ermakov invariant. 



'^-\ 



{i]a — 7]a) 



2 , fV''^ 
a 



co....i(i^)\ (>., 



{f) = 


h 

2m 
hm 

2 




2m 


hm 
2 


{XX*), 


{[x,p\ + ) = 


{xp 


+ px) = 


haa - 


-^^(XX 
- 2di^^^ 



that connects the classical position and momentum (or velocity) and their uncertainties. This 
becomes even more obvious if one realizes that, with the help of a, ip and A and their time- 
derivatives (denoted by overdots), the quantum mechanical uncertainties can now be expressed 
in the following forms 

(20) 

(21) 

(22) 

where x = x — (x), p = p — (p). 

In order to compare the time-dependent kernel ^ with the kernel K{x,x') in ([5]), one must 
take into account that K{x,x') has been obtained via equations ^ and (|3D which describe the 
transformation of x and p into the initial values x' and p' , whereas, K{x, x' , t, t') in ([8]) describes 
the inverse transformation from x' to x. This is expressed, e.g., by the different signs in the 
exponents of ^ and ([8]). In order to make a direct comparison, one must therefore take the 
inverse transformation of ([8]), obtained by changing the sign in the exponent, and interchan- 
ging z and u. Inserting this kernel into equations ([3]) and (JH), one obtains the corresponding 
equations for the time-dependent problem, 

zx — z — = — , (23) 

m ao 

—ux + u — = (24) 

m m 

or, in matrix notation, 

^ -^ ^ ^ ^ ^ mi I ]. (25) 




-u u 



p ^"- \ p 



Again, the transformation matrix M has the determinant zu — uz = 1, which corresponds, 
according to (J17p . to a kind of conservation of angular momentum (in a complex plane). However, 
different from the time-independent case, the initial state is not only characterized by the initial 
position x' and momentum p' but, also, by the corresponding initial uncertainties since oq = 
(^(x^)o)^'^ is proportional to the initial position uncertainty and (for a minimum uncertainty 
wave packet with (x)q(p^)o = ^^/4) the inverse ^ = (;^(p^)o)^ is proportional to the initial 
momentum uncertainty, i.e., ^- oc ,^„, , ^^ oc /„, . 



Wigner Function and Canonical Transformations in Time-Dependent Quantum Mechanics 7 



Following the procedure for the inverted propagator ([8]) outlined in [7], by applying the 
Wigner transformation ([B]), one arrives at a similar result. The kernel K{x,x' ,p,p' ,t,t') (being 
time-dependent via z{t) and u{t)), which provides the Wigner function W{x,p,t) via 



/ + 00 
dx'dp K{x,x' ,p,p ,t,t')W{x' ,p ,t' = 0) 
-oo 



with 



J2 



W{x',p,t' = 0) 



exp 



p. 



7rh I h 



2(x2)o 2(p2)o 

2 



m 



X 

— I + 

"0 



aop 



m 



(26) 



(see, e.g., [H]), is given by the product of two delta functions where now, however, x' is re- 
placed by ^ and p' is replace 
determined by ^, ([SID, i.e., 



placed by — and p' is replaced by ^^^^ and the transformed variables in the delta functions are 



K{x,x' ,p,p ,t,t' = 0) 



X 



zx — z- 



V 



m 



aop 



m 



. P A 

u ux 

m 



Applying this kernel to the initial Wigner distribution function (|26p yields the function 
W{x,p,t) as 



W{x,p,t) = ^exp<^ — - 



zx — z — 1 + I u ux 

mJ V m 



Using the definitions (J13p and (J16p of A and its relation to the uncertainties as given in 
equations (|20p - (j22p . finally allows one to write 



W{x 



1 f 2 r, <- 



)x — {[x,p\-[.)xp + {x )p 



(27) 



where the time-dependence of the uncertainties is determined totally by the time-dependence 
of z{t) and u{t). In the case of time-dependent Gaussian wave packets, the classical time- 
dependence is expressed by the fact that the maximum of the wave packet follows the classical 
trajectory. This is taken into account by shifting the variables of position and momentum 
from X to X = X — (x) and p to p = p — {p). Since (x) and {p) are purely time-dependent 
quantities, x and p can replace x and p in equations ([3]), (jl]) since these equations only contain 
derivatives with respect to space, not time. So, x and p in (f27P would be replaced by x and p 
which would lead to the result obtained in |16j showing the connection between the exponent 
of the time-dependent Wigner function and the dynamical Ermakov invariant that is connected 
with the parameters z and u of the time-dependent kernel K{x, x', t, t') and has been defined in 
equation (fTOl) (for details see also [TT]). 

In the quantum mechanical phase space picture according to Wigner, this results not only in 
changing initial position- and momentum-uncertainties into their values at time t but, also, an 
additional contribution occurs from the time-change of (x^), or a^, respectively, expressed by 
the term proportional to ([x,p]+), or da, respectively, in the exponent of W{x,p,t). 

All these quantum dynamical aspects are contained in the time-dependent parameters z 
and "u, entering the transformation matrix in (j25p . In particular, the change of the position 
uncertainty (proportional to a) is taken into account by the parameter u, which can be expressed 

as [m 



u 



za — zaa 



/ 771 



aoPo 



[f]a — r]aa\. 



(28) 
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For constant uncertainty a = oq, u is simply proportional to the classical velocity r]{t), for 
d 7^ 0, however, the situation can become quite different. As an example the free motion shall be 
discussed briefly, since there a = a{t), which is expressed in the spreading of the corresponding 
wave packet solution. For this purpose, also u is now given in terms of rj and rj where the 
equations of motion (fTSj) (for ??(i)) and (fT8|) (for a{t)) are applied. So it follows from ([28|) that 



/ m 

u= 

\aoPo 



2 1 
rjaa — ri { a -\ ^ 



For constant q = oq, all terms proportional to a vanish and the transformation matrix in ([25 
can be written as 

M=i K -M = r^U? "."V 



v"opoy 



u u J \aoPoJ \ -7 OqT] 



^0 



For the harmonic oscillator with constant width (note: there is another solution with oscil- 
lating width) a = ao = -j= and r]{t) = ^sinwi (??(0) = 0), ?) = VQCOsujt (77(0) = vq = ^), 
M turns into 



( rn \ ( Vo cos cut -^ sin uit \ f -^^ cos ujt -QosinL<jt \ 
\ otQfQ J \ vq sin Lot ^ cos ujt j \ ^ sin ujt uq cos ujt I 

i.e., (up to the constant ao that also occurs in the column vectors) just the classical result is 
reproduced. 

However, for the free motion with 77(f) = vot, 77(f) = vq, for constant a = a^, one would 
obtain 

^ _ ( m \ f vo -Vot \ 
^''~ \aopo) y ^vot alvo J ' 

that is different from the classical situation where the matrix element in the first column and 
second row would be zero. The consequence for a free motion wave packet with constant width 
would be that the transformation matrix would no longer describe a canonical transformation, 
since its determinant would no longer be equal to 1 but 

det(Mfi.) = 1 + ' 

which just describes the time-dependence of the wave packet spreading. 

For d 7^ 0, the well-known time-dependence of the wave packet width, given by a^(f) 
2I 

and obtained as a solution of equation (jlSp for u; = 0, leads to the correct 



2 

"0 



1 + 1:^ 



^0, 
transformation matrix 

\aoPoJ V «o^^o J 

with det(Mfi.) = 1. This shows explicitly the influence of the time-dependence of the uncer- 
tainty a on the transformation describing the dynamics of the system. 

It should also be mentioned that the determinant of M, written in terms of rj, 77, a and a 
takes just the form of the Ermakov invariant, i.e., 

M= ( ^^) ( .. ^^,.2 IN . o"'' ■ I (30) 

" " -?7aa -|- r] (q -|- ^ j 77a — r]aa J 
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yields 



det(M) 



f m \ 
\aoPoJ 

i—Y 






{f]a - ar]Y + ( — ] 



(31) 



At this point it appears appropriate to discuss some differences characteristic of our method in 
comparison with other approaches that try to find a quantum analogue to the classical canonical 
transformations in the case of time-dependent quantum systems, particularly, when the Wigner 
phase space formulation is applied. 

So in [T71 [18] the Wigner function for the free motion and Gaussian-shaped distribution 
function is given by {h = 1) 



W{x,p,t) = — exp 

vr 



p_ 
m 



t] +bp' 



h^ 



with constant parameter b (that would, in our notation, correspond to 6 = ■^(x^)o = 3i ~ 2("^) 
(for a minimum uncertainty wave packet)). The equivalence to the conserved volume element 
in classical phase space is intended to be established by the statement: "this distribution is 
concentrated within the region where the exponent is less than 1 in magnitude" and, for the 
choice of a coordinate system where 6 = 1, "the above phase space distribution function is a circle 
at t = 0. As time progresses, the circle becomes a tilted ellipse while preserving its area." This 
elliptic deformation is a canonical transformation that corresponds to our transformation (j29p 
but, in the case of these authors' description, the time-dependence in their transformation 
only originates from the classical dynamics, transforming the initial position x into x — —t 
(with constant p). The correct spreading Gaussian wave packet for the free motion, however, 
obviously has a time-dependent width, corresponding to a time-dependent parameter b in the 
notation of the references quoted. This still holds if one considers the corresponding Wigner 
function. However, if the width in position space is time-dependent, in our notation a ^ 0, 
a third term in the exponent of the Wigner function, namely, the one taking into account the 
position-momentum-correlations (see equations (j22p and (j27p ). must occur. This term is missing 
in the above-mentioned approach. 

In another approach [19], a relation between the Wigner function for two-photon coherent 
states and the one for Glauber coherent states has been established that looks quite similar to 
our transformations (fTj) or (j25p . However, this method is based on combinations of creation 
and annihilation operators. It has been shown recently [20] that the Ermakov invariant (that 
is, up to a constant factor, equivalent to the determinant of our transformation matrix (j25p : see 
also pip ), can be factorized into two terms that are a kind of (complex) generalization of the 
creation and annihilation operators and turn into these for a = (for further details, see |20|). 
Therefore, also in the approach [19] based on the usual creation and annihilation operators, the 
time-dependence of the uncertainties expressed by d 7^ is not taken into account. 

It is also known in the literature that quantum uncertainties are related with classical error 
margins and that, particularly for quadratic Hamiltonians, the classical error margins satisfy the 
classical Hamiltonian equations [21]. Furthermore, in our case, it can even be shown that the 
quantum uncertainties obey equations of motion that can be derived from a quantum-uncertainty 
Hamiltonian function that now provides the time-evolution of the quantum uncertainties in 
a canonical formalism. This Hamiltonian is nothing but the ground state energy or, respectively, 
the energy contribution of the momentum and position fluctuations to the overall energy of the 
Gaussian wave packet solution for the system. To make these formal aspects obvious, we write 
the energy of the system, calculated as the mean value of the Hamiltonian operator using the 
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Gaussian wave packet solutions of the corresponding time-dependent Schrodinger equation, in 
the form (for the harmonic oscillator) 



{H.p) 



h'^'-' -'^^'<:''-> 






i(p>^ + ^.V->'M + 



- (a +a if ) + -u; a 



i^<^'^-'^^"<''-^ 



In order to establish a Lagrangian/Hamiltonian formalism for the uncertainties, we assume 
that a corresponding Lagrangian C can be written as the difference between kinetic and potential 
energy fluctuations, but now expressed in terms of the variables a, (p and the corresponding 
velocities a, if, i.e. 



C{a, d,ip,(p) = f — V = — (^a + a (f — uj a ) . 
The corresponding Euler-Lagrange equations are then 



d_dC_dC_ ^^_^-n 

dt dip dtp ' dt da da 

From the first equation follows ^7(0'^^'/') = 0' o'^' '^^V' = const, in agreement with equation ()17p : 



from the second equation follows a + J^a 
const = 1). 



(p^a 



3~) 



equivalent to equation (fTHj) (for 



The corresponding canonical momenta are then given by 



dL h 



dC h 



a^ = r^=^-' d^ = r=p- 



With the help of these definitions, the quantum energy contribution E 
written in a Hamiltonian form as 



T + V can be 



n{a,pa,ip,p^) 



p\ 



h ha^ 



h 



Pip ,"'22 

4 



It is straightforward to show [13] that the corresponding Hamiltonian equations of motion 
again reproduce the results (J17p and (jlSp . only now expressed with the help of the canonical 
momenta. An interesting point is that because of relation (J17p . i.e. if = ^i, the canonical 
"angular momentum" p^p has the constant value Pi^ = |, a value that does not usually describe 
an orbital angular momentum but the non-classical angular momentum-type quantity spin. 

Finally, using these results, the uncertainty product can be expressed as 



U = {x'){j}')=p' + {ap^Y. 



In this context it should be mentioned that some authors (e.g., |22 1 I23 ^ [M]) assume that the 
role of the phase space volume in quantum mechanics is played by the square root of the so-called 
"invariant uncertainty product" {Sp'){j?) — j([x,p]+)^ = ^, which is definitely a constant of 
motion since it is equivalent to our quantity p^. So, this square root would just be p^p = i^oP''^ 
which, due to the equivalence between equations Q and (|17p . leads back to the requirement for 
the determinant of our transformation matrix (|25p to be equal to 1. 
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4 Conclusions and perspectives 

In classical mechanics, canonical transformations can be characterized by certain properties, like 
the conservation of a volume element in phase space under these transformations, corresponding 
to Liouville's theorem that can easily be checked by calculating the Wronskian determinant of 
the transformation that should be equal to 1. 

In quantum mechanics, this conservation law appears to have its correspondence in the conser- 
vation of a so-called "invariant uncertainty product" that holds for systems with time-dependent 
and time-independent quantum uncertainties since any explicit time-dependence of the uncer- 
tainties is compensated for by subtracting a term proportional to ([x,p] + )^. The remaining 
conserved quantity corresponds to the conservation of p^p^ the "angular momentum" for the 
motion of A(t) in the complex plane. Due to the equivalence of equations ([9]) and (J17p . this 
conservation law is identical with the requirement that the determinant of our transformation 
matrix (p5|) for the time-dependent quantum problem must be equal to 1. From equations ([30|) 
and (I3ip . it finally follows that this requirement is identical with the existence of a dynamical 
invariant for the system, the so-called Ermakov invariant. 

Another major result of our analysis is that, in the time-dependent quantum mechanical prob- 
lem, the transformation ()25p corresponding to the classical linear canonical transformation ([1]) 
and its time-independent quantum mechanical analogue ([3]), (jl]) does not only transform the 
initial position and momentum into its values at a later time but, also, does the same sim,ul- 
taneously with the corresponding uncertainties! In how far this is connected with the existence 
of a Lagrangian/Hamiltonian formulation of the dynamics of the quantum uncertainties will be 
further investigated. 

So far, the discussion of the time-dependent case included only systems where the potential 
is at most quadratic in its variables. This might not be as restrictive as it seems at first sight 
since one may sometimes perform canonical transformations to reduce a given Hamiltonian to 
a quadratic form [25] which has been shown explicitly by Sarlet for some polynomial Hamil- 
tonians. To what extent this method can also be applied in our case requires more detailed 
studies. 
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